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Abstract 

Anisotropic homogeneous metamaterials that are neither wholly dissipative nor wholly active at a 
specific frequency are permitted by classical electromagnetic theory. Well-established homogenization 
formalisms indicate that such a metamaterial may be realized quite simply as a random mixture of 
electrically small (possibly nanoscale) spheroidal particles of at least two different isotropic dielectric 
materials, one of which must be dissipative but the other active. The dielectric properties of this meta¬ 
material are influenced by the volume fraction, spatial distribution, particle shape and size, and the 
relative permittivities of the component materials. Similar metamaterials with more complicated lin¬ 
ear as well as nonlinear constitutive properties are possible. Dynamic control of the active component 
material, for example via stimulated Raman scattering, affords dynamic control of the metamaterial. 


1 Introduction 

Causality mandates that electromagnetic fields must attenuate as they propagate inside homogeneous, pas¬ 
sive linear materials [1]. Weak electromagnetic fields in certain spectral regimes can be amplified in some 
homogeneous linear materials, provided that strong electromagnetic fields can pump in the energy needed 
for amplification [2, 3, 4, 5]. Conditions on the constitutive parameters of homogeneous linear materials 
have been derived to determine if a linear material is either dissipative or active, but not both, at a specific 
frequency [6]. 

Active component materials feature prominently in the field of metamaterials [7], in order to overcome 
losses [8] and to enhance performance [9]. While anisotropic structures containing active components have 
been reported upon previously [10, 11, 12], the prospect of simultaneous attenuation and amplification of 
electromagnetic fields in homogeneous materials at a specific frequency, depending upon the orientation of 
electric fields, has not been considered hitherto. Materials supporting amplification and attenuation at the 
same frequency are necessarily anisotropic. A host of applications for them can be envisaged. These materials 
could be used in directional coupling devices [13] as well as in spatially discriminatory and/or frequency- 
discriminatory optical amplifiers [14]. Amplification could be dynamically controlled by exploiting, for 
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instance, the phenomenon of stimulated Raman scattering [2, 15], thereby affording dynamic control. That 
is, if the active component material were Raman active, then the degree of amplification achieved for a probe 
laser beam at a desired frequency could be dynamically controlled by means of a strong pump laser beam at 
a determined frequency which induces Raman transitions. The difference between the probe and pump laser 
frequencies is specified by the energy levels of the Raman-active material. The ability to suppress radiation 
leakage in certain directions while promoting propagation in other directions could be harnessed to amplify 
surface-plasmon polaritons and reduce optical noise in biosensing applications [16], for example. Radomes 
for enhancing or reducing directionality of radiation from optical antennas could be made of these materials, 
leading to enhancements in the efficiency of photodetection, light emission, and sensing [17]. 

Motivated by these potential applications, here we propose dynamically controllable anisotropic materials 
which simultaneously exhibit both dissipation and amplification at a specific frequency, depending upon the 
orientation of electric fields. 


2 Simultaneous attenuation and amplification 


Consider a generally anisotropic, homogeneous, linear, dielectric material characterized by the frequency- 
domain constitutive relations 

D(r,cd) = eog{uj) • E(r,cd) 

B(r,cd) = /ioH(r,cj) 

where e{uj) is the relative permittivity dyadic [18, 19] at angular frequency cj, and Cq and //q are, respectively, 
the permittivity and permeability of free space. The time-averaged dissipated power per unit volume is given 
by [19] 

(3(r,w) =-^^E*(r,w) • [|(w) - |(w)] •E(r,w), (2) 

where i = \/— 1 ; an exp(—icdt) dependence on time t is implicit; the superscript * denotes the complex 
conjugate; and e{uj) is the hermitian conjugate of €{uj). 

At a specific angular frequency cd, the chosen material is classified as [ 20 ]: 

• dissipative if Q{r,uj) > 0 , which requires the dyadic i [ 6 ( 61 ;) — e(co’)] to be negative definite; or 



• active if < 0 , which requires the dyadic i \g{uj) — e(c(;)] to be positive definite. 

All eigenvalues of a negative/positive definite dyadic are negative/positive [ 21 ]. Henceforth, for compact 
representation, the dependency of Q and E on cj and r, will not be explicitly stated; similarly, the dependency 
of e (and its components) on uj will not be explicitly stated. 

The foregoing classification fails to accommodate to the prospect of i (e — e) being indefinite [21]—i.e., 
when some but not all eigenvalues of i (e — e) are positive, the remaining eigenvalues being negative. For 
orthorhombic materials [ 22 ], i (e — e) is indefinite provided that Im {e} is indefinite, where the operator 
Im{-} delivers the imaginary part. 

The simplest material, at least from a mathematical perspective, for which Im {e} is indefinite is a 
uniaxial dielectric material whose relative permittivity dyadic has the form [19] 

e = (/ — uu) + el'uu, (3) 


with Im Im {ell} < 0 . Herein / is the identity dyadic and the unit vector u is 
optic axis. For this material, it may be inferred from Eq. (2) that Q for E directed 
sign to Q for E directed perpendicular to u, because 


Q = ^. 
^ 2 L 


Imje-^jE 


(/ — uu) • E 


-fim 


E* • uu • E 


parallel to the material’s 
along u has the opposite 


( 4 ) 


Thus, there is dissipation associated with certain orientations of the electric field but amplification with 
other orientations. 
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3 Realization as a homogenized composite material 

Is it possible to realize a uniaxial dielectric material for which Im {e} is indefinite? We now demonstrate, 
using well-established theoretical formalisms based on the homogenization of particulate composite materials, 
that materials with indefinite Im {e} may be conceptualized as homogenized composite materials (HCMs). 
As these engineered materials will simultaneously exhibit both amplification and dissipation at the same 
frequency, they should more properly be called as metamaterials [23]. In the visible spectrum (380-770 
THz), the maximum linear dimensions of the component particles are required to be less than ^30 nm [24]. 

Consider a composite material that is a mixture of two component materials labeled ‘a’ and ‘b’. Dispersed 
randomly as identically oriented, electrically small, conformal, spheroids, both component materials are 
isotropic dielectric materials with relative permittivities and e^. The surface of a spheroid, relative to its 
centroid, is prescribed by the position vector r = • r^, where p > 0 is a linear measure of particle size, 

the shape dyadic 

U. = —^ (I - uu) + U uu, {U > 0) (5) 

yU ~ 

contains the shape parameter I/, and the unit vector prescribes the surface of the concentric unit sphere. 
Oblate spheroidal particles are characterized by U G (0,1), prolate spheroidal particles hy U > 1, and 
spherical particles hy U = 1. The volume fraction of component material ‘a’ is fa, while that of component 
material ‘b’ is /^ = 1 — /q. The spheroidal shape and the identical alignment of the component particles 
endows the HCM with uniaxial symmetry [18, 19]; i.e., the relative permittivity dyadic of the HCM has the 
form 

Ihcm = ^hcm (I - u u) + 4cmU u- (6) 

Porous columnar-thin-film sections and nanoparticle arrays can be fabricated to realize such HCMs [25, 26]. 

4 Estimates of constitutive parameters 

Let US now present representative numerical estimates of as yielded by the Bruggeman formalism 

[27], the strong-permittivity-fluctuation theory (SPFT) [28, 29], and the Maxwell Garnett formalism [30, 31]; 
these estimates are identified by replacing the subscripts ‘HCM’ by ‘Br’, ‘SPFT’, and ‘MG’, respectively. 
For the purpose of illustration, suppose that the component material ‘a’ is an active material with 
= 2 — 0.03L This value of lies comfortably within the range typically employed for active components 
of metamaterials in the visible regime. For example, a mixture of two commonly used amplification materials, 
namely Rhodamine 800 and Rhodamine 6G, possesses a relative permittivity with imaginary part in the range 
(—0.15, —0.02) and real part in the range (1.8, 2.3) across the frequency range 440-500 THz, depending upon 
the relative concentrations and the external pumping rate [11]. Component material ‘b’ is taken to be a 
dissipative material specified by = 3 + 0.05i, and the shape parameter U = b. The real and imaginary 
parts of the HCM’s relative-permittivity scalars and as estimated using the Bruggeman formalism 
[27], are plotted against the volume fraction fa in Fig. 1. The real parts of and decrease in an 
approximately linear manner from 3 to 2 as fa increases from 0 to 1. Both Im and Im decrease 

uniformly from Im to Im {ea} as fa increases from 0 to 1, but the decrease in Im is more distinctly 

nonlinear than that in Im | |. The data show that 

(i) for 0 < fa < 0.52, the HCM exhibits only dissipation, regardless of the orientation of the electric field 

E; 

(ii) for 0.63 < /a < 1, the HCM exhibits only amplification, regardless of the orientation of E; and 

(iii) for 0.52 < fa < 0.63, the HCM exhibits simultaneously both dissipation and amplification, depending 
on the orientation of E. 
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Bruggeman formalism 




Figure 1: (Color online) Real and imaginary parts of (red, solid curve) and (blue, dashed curve) 
plotted against the volume fraction /a, when = 2 — 0.03i, 65 = 3 + 0.05i, and U = b. The yellow rectangle 
indicates the /o-range where Im < 0 but Im < i > 0 . 


Parenthetically, the HCM represented in Fig. 1 is not an ‘indefinite’ material [18] since Re is positive 

definite. 

The shape of the component particles plays an important role in the simultaneous exhibition of dissipation 
and amplification. In Fig. 2, the real and imaginary parts of and are plotted against 1/, when 
fa = 0.535 and the relative permittivities of the component materials are the same as for Fig. 1. The 
component particles are: oblate spheroids for 1 / < 1 , prolate spheroids for 1 / > 1 , and spheres for 1 / = 1 . 
Both the real and imaginary parts of decrease nonlinearly as U increases whereas both the real and 
imaginary parts of increase nonlinearly as U increases, with at = 1. In the vicinity of 

= 1 , Im{e;^^} and Imje^^j are both positive. However, for U < 0.75, we find that Imje^^j < 0 but 

Im{e;g^} > 0. Incidentally, the Bruggeman formalism for a particulate composite material yields the same 
results in the limit U ^ 0 as for a periodically laminated composite material [32]. Also, for U > 2.4, we 

find that Imje'^^j > 0 but Im{e;g^} < 0. Thus, the HCM exhibits simultaneously both dissipation and 
amplification provided that the component spheroidal particles are either sufficiently flattened or sufficiently 
elongated. 

The effect of the relative permittivity of the active component material on the simultaneous exhibition 
of dissipation and amplification is taken up through Fig. 3. Therein, the real and imaginary parts and 

as well as the product Im Im are plotted against the real and (negative) imaginary parts 

of ea, the product being negative if and only if Im and Im have opposite signs. For calculating 

these results, we fixed = 3 + 0.05i, fa = 0.25, and U = b. The ea-regime for which the HCM exhibits 
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Figure 2 : (Color online) Real and imaginary parts of (red, solid curve) and (blue, dashed curve) 
plotted against the shape parameter [/, when = 2 — 0.03i, 65 = 3 + 0.05i, and fa = 0.535. The two yellow 
rectangles indicate the [/-ranges where the product Im | < 0 - 


simultaneously both amplification and attenuation is characterized by relatively small values of Re{ea} 
and relatively large values of —Im{eo}. Indeed, as —Im{ea} approaches 0, the HCM is either exclusively 
dissipative or exclusively active, this state being attained at larger values of — Im {eo} as Re { 6 ^} approaches 

2 . 

Next let us turn to a more sophisticated homogenization formalism: the SPFT [33]. Unlike the Bruggeman 
formalism, the SPFT can accommodate a comprehensive description of the distributional statistics of the 
component materials, via the characteristic functions 


^e{r) 


1 , 

0 , 


rev, 


r^Ve, 


i e {a,b}. 


(7) 


Herein, the regions occupied by the component materials ‘a’ and ‘b’ are identified by V and V, respectively. 
The ensemble average ($^(r) )e equals the volume fraction fe, i e {a, 6}. In the usual implementation of the 
bilocally approximated SPFT [33], the distributional statistics of the component materials are characterized 
by the second moment 




/£, 1^-' • (r-rO| <F, 
/I, 1^-' • (r-rO| >F, 


(8) 


for I G {a, 5}, with L as the correlation length. 

The real and imaginary parts of the HCM’s relative-permittivity scalars and as estimated 

using the bilocally approximated SPFT, are plotted against the normalized correlation length k^L in Fig. 4, 
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plotted against the real and (negative) imaginary parts of the relative permittivity when = 3 + 0.05i, 
fa = 0.25, and [7 = 5. 


where ko = For these calculations, we set = 2 —0.03i, = 3 + 0.05i, fa = 0.535, and 7/ = 5; also, 

we used an extended version of the SPFT [34] which explicitly accommodates the particle-size parameter p. 
Results are presented in Fig. 4 for p G {0, 0.51/, L}. In this figure, the real and imaginary parts of both 
and e|ppp increase uniformly as k^L increases, for all values of pjL considered. Also, Im |e|ppp| > 0 for 

all values of koL and p/L. However, Ini{e; 5 ppp} < 0 for low values of koL while Im{e; 5 ppp} > 0 for 
high values of koL. Furthermore, the transition from negative Imje^ppp} to positive Imje^ppp} occurs 
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at lower values of k^L when the size parameter p is larger. Thus, the HCM exhibits simultaneously both 
amplification and dissipation provided that both L and p are sufficiently small. 


SPFT 



/CqL 

SPFT 



/cqL 

Figure 4: (Color online) Real and imaginary parts of (thick curves) and e|ppp (thin curves) plotted 

against the normalized correlation length koL^ when = 2 — 0.03i, = 3 + 0.05i, fa = 0.535, and U = 5. 

The size parameter p = 0 (red, solid curves), 0.5L (blue, dashed curves), and L (green, broken dashed 
curves). 

Lastly, we present estimates of the HCM’s relative permittivity dyadic provided by the Maxwell Garnett 
formalism [27], whose provenance is quite different from that of either the Bruggeman formalism or the 
SPFT. We consider that particles of component material ‘a’ are dispersed randomly in the component 
material ‘b’. The Maxwell Garnett formalism is restricted to dilute composite materials {fa ~ 0.3). The 

real and imaginary parts ej^Q and as well as the product Im Im are plotted against 

the real and (negative) imaginary parts of in Fig. 5. The component materials are as specified for Fig. 3, 
with one exception: the shape parameter U applies only to component material ‘a’ since particle shape is 
irrelevant to the host material in the Maxwell Garnett formalism. The extent of the e^-regime for which 
Im Im < 0 in Fig. 5 is slightly smaller that the analogous ea-regime in Fig. 3. But, other 

than this relatively minor difference, the plots of Im {ej^o} Im in Fig. 5 and of Im Im 

in Fig. 3 are very similar. 

The estimates of f provided by the Bruggeman and Maxwell Garnett formalisms—as represented 
in Figs. 1, 2, 3, and 5—are in close agreement over the volume-fraction range appropriate to the Maxwell 
Garnett formalism (i.e., 0 < ^ 0.3). The estimates provided by the two formalisms are identical in 

the limit /a ^ 0 and very small differences emerge as fa increases. Indeed, the corresponding plots of 
the real and imaginary parts of the components of and in Figs. 3 and 5, respectively, are almost 
indistinguishable to the naked eye. These estimates are also in close agreement with the corresponding 
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Figure 5: (Color online) Real and imaginary parts ej^Q and as well as the product Im Im 

plotted against the real and (negative) imaginary parts of the relative permittivity e^, when = 3 + 0.05i, 
fa = 0.25, and U = b. Particles of component material ‘a’ are taken to be dispersed randomly in the 
component material ‘b’, and the shape parameter U applies only to component material ‘a’. 


estimates provided by the bilocally approximated SPFT, as represented in Fig. 4. The imaginary parts 







of the components of f estimated by the SPFT deviate slightly from their corresponding Bruggeman 
and Maxwell Garnett counterparts, with the deviation increasing in magnitude as the correlation length L 
and the size parameter p increase. This deviation reflects the fact that the SPFT formalism accommodates 
coherent scattering losses via L and p [33], whereas the Bruggeman and Maxwell Garnett formalisms do not. 

Rigorous theoretical bases have been firmly established for each of the three homogenization formalisms 
employed here [27, 28, 29, 30, 31, 33]. Techniques based on such homogenization formalisms have the 
advantages over full-wave numerical techniques, based on the finite-element method or the finite-difference 
time-domain method [35, 36, 37], for example, that they provide estimates of the constitutive parameters 
which are independent of the shape and size of the bulk material involved, and these estimates apply for 
all possible incident fields with sources not located in the bulk material. However, it is important to bear 
in mind that the predictions of constitutive parameters provided by any homogenization formalism are 
estimates. The ultimate checks on such estimates can only be provided by careful experimental studies. 


5 SIMULTANEOUS ATTENUATION AND AMPLIFICATION 
Exemplified 

Power flow associated with electromagnetic propagation is represented by the time-averaged Poynting vector 
P = (1/2) Re {E X H*}, with the operator Rej-} delivering the real part. Since V • P = — Q in a region 
devoid of externally impressed sources [1, 6], amplification and dissipation should be discernible through 
plane-wave propagation. 

Suppose that the direction of propagation is parallel to the unit vector a. An ordinary plane wave will 
propagate in the chosen HGM with wavenumber kor = ^hcm^ regardless of the angle 0 = cos“^ (a • u). 
An extraordinary plane wave will propagate in the chosen HGM with wavenumber 


^ex 



e 


± e" 
HCM^HCM 


(9) 


that does depend on 0. 

Gonsider two examples: (i) When fa = 0.6 in Fig. 1, = 2.353 — 0.006i and = 2.398 -t- 0 .002i, 

which yield kor = (1.534 — 0.002i) ko with Ke{kex} > 0 for all 0 and lm{kex} > 0 for 60.5° < 0 < 119.5°. 
Hence, the ordinary plane wave is amplified for all propagation directions while the extraordinary plane wave 
is amplified for 0 G (0°,60.5°) U (119.5°, 180°) but attenuated for 0 G (60.5°, 119.5°). (ii) When U = 0.25 in 
Fig. 2, = 2.456-F0.006i and = 2.383 — 0.006i, which yield kor = (1.567 + 0 .002i) ko with Re {kex} > 0 

for all 0 and lm{kex} < 0 for 44.1° < 0 < 135.9°. Hence, the ordinary plane wave is attenuated for all 
propagation directions while the extraordinary plane wave is amplified for 0 G (44.1°, 135.9°) but attenuated 
for 0 G (0°, 44.1°) U (135.9°, 180°). 


6 Closing remarks 

In conclusion, according to the estimates afforded by three different well-established homogenization for¬ 
malisms, an HGM that exhibits simultaneously attenuation and amplification of electromagnetic fields at a 
specific frequency may be realized quite simply as a random mixture of electrically small spheroids of two 
different materials. Both component materials are isotropic dielectric materials, one of which is dissipa¬ 
tive while the other is active. The realization of such an HGM depends upon the volume fraction, spatial 
distribution, particle shape and size, and the relative permittivities of the component materials. Dynamic 
control of the active component material, for example via stimulated Raman scattering, affords dynamical 
control of the HGM. Thus, a new class of metamaterials is proposed which are neither wholly dissipative 
nor wholly active. Although we have illustrated the concept with uniaxial dielectric HGMs, particulate 
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composite materials with more complicated linear and/or nonlinear constitutive properties and displaying 
both dissipation and amplification at the same frequency can be designed [18]. Finally, even periodically 
laminated composite materials [32] offer similar promise, but particulate composite materials may be more 
readily fabricated than their laminar counterparts. 
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